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Abstract 

The solutions of the time independent Schrodinger equation for non- 
Hermitian (NH) Hamiltonians have been extensively studied and calcu- 
lated in many different fields of physics by using L 2 methods that origi- 
nally have been developed for the calculations of bound states. The exist- 
ing non-Hermitian formalism breaks down when dealing with wavepack- 
ets(WP). An open question is how time dependent expectation values can 
be calculated when the Hamiltonian is NH ? Using the F-product for- 
malism, which was recently proposed, [J. Phys. Chem., 107, 7181 (2003)] 
we calculate the time dependent expectation values of different observable 
quantities for a simple well known study test case model Hamiltonian. We 
carry out a comparison between these results with those obtained from 
conventional(i.e., Hermitian) quantum mechanics (QM) calculations. The 
remarkable agreement between these results emphasizes the fact that in 
the NH-QM, unlike standard QM, there is no need to split the entire space 
into two regions; i.e., the interaction region and its surrounding. Our re- 
sults open a door for new type of WP propagation calculations within 
the NH-QM formalism that until now were impossible. In particular our 
work is relevant to the many different fields in physics and chemistry where 
complex absorbing potentials are introduced in order to reduce the prop- 
agation calculations into a restricted region in space where the artificial 
reflections from the edge of the numerical grid/box are avoided. 



1 Introduction 



In the last two decades there has been an increasing interest in non-Hermitian(NH) 
quantum mechanics (QM) in many different fields of physics. In many stud- 
ied cases in atomic, molecular and nuclear physics the Hamiltonian is NH 
due to the use of different type of analytical continuation transformations. 
One such method is complex scaling (CS) where H(r) -> H(re l9 ) Hj Ej |H • 
Another common method is the introduction of complex absorbing local en- 
ergy independent potentials (CAPs), also known as optical potentials, into the 
Hamiltonian |5] [5] such that H — > H + Vcap- The complex eigenvalues and 
eigenfunctions of the NH Hamiltonians, in these studies, were calculated using 
square-integrable methods that were originally developed for the calculation of 
bound states (see references to reviews mentioned above). The need to replace 
the scalar product which is used in conventional QM by a generalized inner 
product is well understood and has been discussed quite extensively in the lit- 
erature. See for example the different approaches for generalized inner products 
in NH-QM in Ref. [7] and Ref.jS] and references therein. However, for the 
sake of clarity in Section 12. II we provide a brief account of the generalized inner 
product (so called c-product 9 ) which has been used until now for calculating 
time independent observables in NH-QM0]. It is quite straightforward to ap- 
ply the same generalized definition of the inner product to NH time periodic 
Hamiltonians ( see for example Ref.0]), which describe the interaction of atoms 
or molecules with oscillating electric field. The application of the NH formalism 
to laser driven system has been very useful in explaining physical phenomena. 
For example, remarkable agreement has been achieved between the results ob- 
tained from theoretical NH-QM calculations and the experimental results of the 
probability of generating odd high-order harmonics (i.e., generation of high en- 
ergy photons resulting from the absorbtion of many low energy photons) 10 . 
However, one important problem remains open. How to calculate the entire 
high-order harmonic generation spectra and not only the probabilities to obtain 
even or odd harmonics ? This question is related to a more general question. 
How to calculate time dependent expectation values in NH-QM ? In NH-QM 
the propagated right wavepacket (WP) is the solution of the time dependent 
Schrodinger equation (TDSE) where H is a NH Hamiltonian. 

&*(r,t)=»ft^*(r,t) . (1) 

Extensive WP propagation calculations were carried out for complex potential 
energy surfaces by solving the NH-TDSE (see for example Refs. |lll 1121 1131 H^l 
115)). However, in order to calculate time dependent expectation values one 
needs to carry out WP propagation calculations from t = to t = -co. The 
problem is that in NH-QM it is impossible to propagate the initial WP from t=0 
to t = -co. We briefly review this known problem in Section 12.21 and relate it 
to the concept of time asymmetry in QM which is a subject of numerous studies 
in the literature (see Refs. ^3 El \IM Ell an d references therein) . 
Very recently a possible solution to the time asymmetry problem in NH-QM 
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was proposed by Moiseyev and Lein[2n|- The solution requires modification of 
the currently used inner product. As we will show in Section[2]this modification 
of the well known generalized inner product (briefly reviewed in Section I2.1fl 
which we will term "F-product" (FP) is applicable for cases where the Hamil- 
tonian is either a time-independent or time dependent NH operator. Naturally 
this modification of the generalized inner product reduces to the conventional 
scalar product when the Hamiltonian is Hermitian. Moiseyev and Lein (2*U| 
have shown that using the modified generalized definition of the inner prod- 
uct not only solves the time asymmetry problem in NH-QM , but also has 
a very important physical implications. They have shown that the analytical 
expression for the high-order harmonic generation spectra (HHGS) which has 
been derived within the framework of the FP modified NH time dependent for- 
malism, clearly indicates why only odd high-order harmonics were observed in 
experiments where high- intensity laser pulses have been used. Moreover, the 
analysis of the expressions they obtained for the HHGS indicates under what 
conditions also even high-order harmonics would be observed. Simulation calcu- 
lations within the framework of the conventional (hermitian) QM based on their 
finding supported their results. The possibility of generating also the even-order 
harmonics as well as the odd ones should help in generating ultra short high 
intensity laser pulses. However, so far the FP generalized inner product as first 
proposed in ref(20| has never been tested for numerical NH calculations. The 
purpose of this work is to calculate time-dependent observable quantities for a 
well known test case study problem within the framework of the FP formalism 
and compare it with the results obtained from the conventional (i.e., Hermitian) 
QM. As we will show here such a comparison clearly shows that time dependent 
observables should be calculated within the framework of the FP formalism and 
also give some physical interpretation of this formalism. 



2 Non-Hermitian Hamiltonians, and their con- 
sequence 

Non-Hermitian Hamiltonians have been used to describe a variety of physical 
phenomena. The methods by which the Hamiltonian becomes NH are also 
diverse. Such methods include for example the addition of complex absorbing 
potentials to the Hamiltonian or the scaling of the coordinate by a complex 
factor. The use of NH Hamiltonians results in complex eigenvalues, 

E a =e a - iT a /2 , (2) 

where often the real part, e Q = Re(E a ), is associated with the energy while 
the complex part is related to the decay rate, T a of a metastable state by: 
r a = —2Im(E a ) such that the lifetime of the metastable state is given by h/T a . 
When using NH Hamiltonians one needs to address the aspects that differ from 
conventional QM. 
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2.1 Brief account of the generalized inner product: the 
c-product 

As a result of the NH nature of the Hamiltonian, there is a need to define a 
generalized inner product. The question of the definition of the inner product 
when the Hamiltonian is NH is a crucial one, especially in time dependent 
calculations. Let us consider a general case where the Hamiltonian, H can be 
either Hermitian or NH Hamiltonian. Since the concept of a generalized inner 
product is well defined in linear algebra for general not necessarily Hermitian 
matrices |21j we begin our brief review by the representation of the Hamiltonian 
by a matrix H. The Hamiltonian matrix elements are defined as 

Hi j =< ipi\H\tpj > (3) 

where {<fi} are orthonormal square integrable basis set functions in the Hilbert 
space. We assume here that we can truncate the number of the basis functions 
and represent the Hamiltonian by a finite general complex matrix. The matrix 
can be as large as one wishes. The matrix H has right and left eigenvectors. Let 
us denote the right (column) eigenvectors by and the left (row) eigenvectors 
by ($t,y. That is, 

H$f = E a $* (4) 

and 

(*£)*H = 25 a ,(*£)* (5) 
By taking the transpose of EqJSJone gets, 

H*$£, = E a ^ L a , (6) 

The matrix H and its transpose H* support the same spectrum. We define the 
generalized inner product, termed c-product [H], such that for non degenerate 
states where E a i ^ E a , 

($ a i|$ a ) = (*£)*• *2 = (7) 

From EqOilit is obvious that when H is Hermitian matrix (i.e., H* = H*) then 
^a' = O^o')* an d the c-product (...|...) is equal to the known scalar product 
< ...|... >. If for example H is a symmetric matrix (note any matrix can be 
transformed to a symmetric form and therefore a complex symmetric matrix is 
a general representation of a NH matrix) then, 

*o = *a (8) 

Do the eigenvectors of a non-hermitian matrix H form a complete set ? A 
complete set implies that the number of linearly independent eigenvectors is 
equal to the dimension of the matrix. It may happen that due to a coales- 
cence of two eigenvectors (coalescence of more than two eigenvectors is very 
unlikely phenomena |2 1 | [2~2 ] ) the number of the linearly independent eigenvec- 
tors is smaller than the dimension of the matrix and the spectrum is incomplete. 
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This coalescence of the eigenvalues happens at a — ab (where "b" stands for a 
branch point in the complex energy plane |2*2*]) for a — > a& the specific eigenvec- 
tor denoted by a = at, is "self-orthogonal" in the sense that, ($ Qb |$ a J = 0. 
However, due to round off numerical errors it is quite impossible to get the 
incomplete spectrum and always in the numerical calculations the eigenvectors 
are normalizable such that, (^I^a) = 1 ESI- As the number of the basis 
functions, {}Pi}, which are used to construct the Hamiltonian matrix in Eq|3J)is 
increased, the function <&^{N) given by: 

N 

<f*(A0=£[$f]^ (9) 

i=l 

gives a better description of the a — th eigenfunction of the Hamiltonian. It 
is out of the scope of this paper to describe the complex (i.e., NH) variational 
theorem that has been proved for complex symmetric Hamiltonian matrices (see 
Ref.0] and references therein). However, on the basis of the complex variational 
principle the exact eigenfunction, cf>^ can be expressed as, 

= lim $f (AO (10) 

N— >oo 



and similarly, 

iv 

N^oo JV->oo- 



lim $£(JV) = lim (11) 

i=l 

where, 

iR tp iR 



H<P« = E a <% (12) 
and 

fr<fe* = E*^* (13) 
Here W is the hermitian conjugate of H such that, 

< V>i \H^\<p j >=[< V , j \H\ t p i >]* (14) 

The basis functions are square integrable functions as described above. Now, 
the generalized inner product (so called c-product) for the left and right eigen- 
function, {4>ai4 > a}i wu l be defined by: 



f (fy&dv = 6 a ,, a (15) 

J all — space 



Note that the c-product, (<f> a \4> a ), may be complex and may have negative real or 
imaginary parts, therefore the c-product is not a metric scalar product [§]. For 
a physical interpretation of the complex density probability <Pa < f ) a see Ref . |24| . 
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2.2 Brief account of the time asymmetry problem in NH- 
QM 

Time asymmetry in physics is a concept closely related to irreversibility which 
has been a subject of many theoretical studies. The conventional Hermitian 
QM is time-symmetric in the sense that it is described by an equation symmet- 
rical with respect to time and by time-symmetric boundary conditions 16]. The 
process of the decay of a metastable state is an irreversible one in the sense that 
it can only be described from a certain time t = to time t — oo. Authors 
such as Bohm [Ijg\, Prigogine |17|. Sudarshan et. al. JH| have constructed for- 
malisms which incorporate the irreversible nature of the resonance phenomena 
into QM usually by relying on the introduction of rigged Hilbert spaces. This 
leads to semi group evolution which distinguished between "prepared" stated 
and "measured" states U^J. However, it has been argued by Nicolaides \Hj\ that 
time asymmetry results from complex energy distribution without the need to 
introduce rigged Hilbert spaces. 

Indeed, in the spirit of Ref.^Hji we will show that in NH-QM, within the spec- 
tral representation of the NH Hamiltonian, time asymmetry will pose a problem 
in the propagation of WP's. The need to carry out WP calculations within the 
framework of NH-QM is a crucial point in the study of systems where the dynam- 
ics are not controlled by a single resonance state. When one tries to propagate 
a WP in NH-QM the evolution of each of the stationary eigenfunctions in Eall2l 
is governed by the NH-TDSE and the stationary solutions of EqQ] are given by: 

i>*(t) = <aqp(-%E a t/h)<jg . (16) 

When one attempts to propagate the corresponding stationary solutions of the 
transposed NH Hamiltonian IP = W* (see Ea llHI) . at time t the left eigenfunc- 
tions will be given by: 

^(i) = exp(+iE a t/K)^ . (17) 

Since E a — £ a — iT a /2 the left eigenfunctions diverge exponentially as t —> oo. 
In Hermitian QM ip^ — i\)* a and E a is real, thus, ip^ does not diverge. Since 
the only difference from the evolution of ip^ is in the + sign in the exponent it 
is similar in a way to propagating the initial state to a negative time. Since in 
NH-QM the energy is complex, the wavefunction cannot be propagated back- 
ward to — oo but only forward to +oo. This time asymmetry problem dissolves 
when studying cases where only one metastable state dominates since when cal- 
culating any observable quantity the time dependent phases of the right and left 
eigenfunctions, exp(±iE a t/fi), cancel out. When one studies the evolution of 
a WP, the right function obeys the NH-TDSE in Eq^and can be constructed 
as a linear combination of the eigenfunctions of the NH Hamiltonian given in 
Ea ll6l The left function will be a linear combination of the eigenfunctions given 
in Ea ll7l Now, when calculating any observable quantity, which does not com- 
mute with H, there will be cross terms between different eigenfunctions with 
time dependent phase of the form, exp(i(E a —E a i)t/K). These terms will diverge 
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at long times as exp((T Q — Y a i)t/2) and prevent the calculation of observable 
quantities of a WP in NH-QM. 

This is, in short, the time asymmetry problem in NH-QM. Due to this problem 
the calculation of observables in NH-QM is feasible only when one long lived 
metastable state controls the dynamics of the system and ipa(t)ipa(t) = <j>a^a 
is time independent. A definition of an observable for a general superposition 
of resonance states is required. This subject will be discussed in Sections |3] and 

El 

3 Wavepacket propagation in the F-Product for- 
malism 

As stated earlier the time asymmetry problem prevents the propagation of WP's 
backward in time in NH-QM. The existing definition of the left eigenfunctions 
of the NH Hamiltonian is based on the c-product as defined in section |2~T1 by 
Ea ll7l which diverges exponentially in time due to the complex part of the 
energy. The one-state probability density ipa^a = ^a&a i s complex and time 
independent. However, as mentioned before the imaginary part of the complex 
energy is associated with the decay of a metastable state therefore one would 
expect that the probability to find the particle in such a state should decay in 
time with a rate of decay r Q . 

In order to impose this decay behavior of a metstable state we apply the finite- 
range F-product (FP) formalism which was proposed in Ref. 20] but was never 
tested by carrying out WP propagation calculations within the framework of 
NH-QM. For the sake of clarity let us first explain the idea behind the F-product 
formalism. As we have seen in section l2~2l when the left function is stationary 
solution of Eq. El h is defined by EqEI Although ipa(t) diverges exponentially 
as time goes to infinity, the norm is preserved and (4> a (t)\tp a (t)) = (<p a \4>a) = 1- 
This result is due to the fact that we integrate over the entire space and the 
probability to find the particle somewhere in space is unity. In order for the 
probability density to decay with time we need to divide the entire space into 
two parts. One part is the interaction region which we associate with our system 
that in time breaks up into sub-systems. The complimentary part is defined as 
a "surrounding" or as the "environment" of our system. We require from an 
isolated single metastable state to decay exponentially (a first order 'reaction') n 
time as the sub-systems escape from the interaction region to the " environment" . 
An acceptable way to separate the "system" from the "surrounding" is by the 
well known Feshbach formalism , which is often used for describing systems in 
nuclear physics. Here the Hamiltonian is split by two operators Q and P which 
project the system into subspaces of discrete and continuum states respectively. 
The resulting Hamiltonian for the bound part of the system is given by: 

HfF — Hqq + VQpGpVpQ , (18) 

where, Hqq is the Hamiltonian for the "system" and Gp = lim £ _» + (E — Hpp + 
ie)^ 1 is the Green operator for the particle in the continuous or dissociative 
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part of the spectrum which describes the " surrounding" to which the particles 
decay. The effective Hamiltonian as defined in Eq. ^|is a NH operator. It is 
clear that here the effective Hamiltonian is NH due to splitting of space into 
an interaction region where our system is located and a "surrounding" which 
absorbs the emitted particles (sub-systems). However, we wish to associate the 
decay phenomena with particles which "escape" from the interaction region to 
its "surrounding" even in cases where the NH formalism takes the entire space 
into account. How to split the entire space in the most general case into an 
interaction region and its "surrounding" ? We will return to this question in 
the next section. Here we argue that when the dynamics is controlled by a 
single isolated metastable state within the NH-QM formalism then we may use 
the finite box quantization approach where the interaction region can be defined 
as a box as large as one wishes. In such a case the integrals are calculated not 
over the entire space, i.e., from — oo to +00, but from —a/2 to +a/2 where the 
size of the box, a, is as large as one wishes, hence the name - finite-range given 
to this new formalism. Since we want to describe a decay phenomena we now 
define the evolution of the left eigenfunctions as, 

ip% = exp(tE£t/7i)<f>a = exp(ie a t/h) exp(-T a t/(2n))^ ■ (19) 

Now the probability do detect the particle in a decaying state will decrease 
exponentially in time according to: 

(V'alV'a) = exp(-T a t/h){4> a \<f> a ) = exp(-T a t/h) . (20) 

This definition also prevents the divergence in time of the left functions and 
enables the propagation of a WP in time. Within the FP formalism a WP will 
be now defined by a superposition of the eigenfunctions of the NH Hamiltonian. 

¥&, = 53c£exp(-iE a i/&)^? ; ^ = 53c^exp(iE*t/n)^ , 

a a 

(21) 

where 

C« = (<f> a \V FP (t = o)) ; c£ = (*fp(* = 0)|&,) ■ (22) 

When the Hamiltonian is NH due to the application of one of the complex scaling 
(CS) similarity transformations, then the initial state ^>Fp(t = 0) is given by 
S^(t = 0) where S is one of the CS transformations . On the other hand 
when CAP is added to the Hamiltonian and the initial state is localized in the 
interaction region it remains unsealed. An important point in the new definition 
of the inner product is that now the right- WP satisfies the NH-TDSE as defined 
in Eq^ whereas the left-WP does not satisfy a corresponding wave equation. 
This is a very important point in our search for solution of the time asymmetric 
problem in NH-QM. The FP definition of the inner product generalizes only the 
time dependence of the c-product definition. Therefore, when the functions are 
time independent the F-product reduces to the c-product. Based on the Eq. 
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the decay of a WP will now be given by: 

(^ F p\^ F p)=^2C^exp(-T a t/h). (23) 

a 

Note that although there are no cross terms between different a's there are still 
interference effects since C^C^ have complex values. Nevertheless the absolute 
value of this expression decays continuously as time passes. The goal of this 
paper is to check the FP formalism for the first time for a simple test model 
problem and compare the results with those obtained from conventional QM 
calculations. Such a comparison is made in section 0] 

4 Application to a simple one dimensional prob- 
lem 

We would like to implement the new formalism on a simple one dimensional 
time independent Hamiltonian, which will serve as a test model: 

^)-^ + (y-0.8)exp(-0.1 a; 2 ) , (24) 

This model Hamiltonian has been often used as a test case for new theories and 
computational methods (see for example H3). The potential in Eq|^]consists of 
a potential between two barriers (See Fig^}. In the example presented below the 
Hamiltonian becomes NH upon complex scaling of the coordinate by a complex 
phase, 9, such that x — > xexp(i9) and is given by: 

^ (x) = e " 2ie £ + (f ~T"" - 8)exp( "°- le2lV) ' (25) 

The model potential shown in FigQ supports a bound state and two metastablc 
resonances below the top of the barrier as well as several other resonances over 
the barrier. The resonances are localized states that decay in time and cannot 
be described by a single eigenfunction of the Hermitian Hamiltonian in Eg 1241 
however in NH-QM they are represented by a single square integrable eigen- 
function of the NH Hamiltonian in Eq[2Hlas can be seen in Fig^ Our objective 
is to show that the FP formalism will yield similar results to those obtained by 
conventional Hermitian QM, without the need to divide our space to "system" 
and "surrounding". For this purpose lets define now the interaction region as 
the area between the top off the two barriers in Fig^ and return to this point 
later on. 

In order to observe the decay of a WP with time we will place a gaussian WP 
in the center of the potential in Eq[^] of the form: 

m = jJw iexp {~i +tkox ) ' (26) 
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where a is the width of the WP and ko its initial momentum. Upon CS the 
right initial state WP becomes: 



* fl (0) = 7^T7iexp(-^+ze^ x) , (27) 



,2)1/4 eX P ^ 2a 2 

while the left initial WP is given by 



pie / lie 2 \ 

^(0) = ^ T7I exp(-^-^%x) . (28) 

Note that at time t = the definitions of the c-product and the F-product are 
identical and (# L (0)|# fl (0)) = 1. The norm of the WP, N F p{t), within the FP 
formalism, as defined in Eg 1211 and is given by: 



FP 



roc 

(*)=/ *%p{x,t)*f P (x,t)dx = J2cZc2exp(-r a t/h) . (29) 



Assuming the decay of the WP is a first order process the effective decay rate 
will be given by: 

k F p{t) = -j t In N FP {t) . (30) 

When the WP populates several resonance states the effective decay rate will be 
time dependent, but when only one resonance state is populated the anticipated 
constant effective decay rate will be obtained, kpp = T res /h. In contrast, the 
decay will not be observed in Hcrmitian QM where the norm is conserved. 
Therefore we choose a new quantity for comparison with the norm in Ea l29l 
which we will label Mqm and is defined by the part of the WP which is localized 
inside the interaction region. 

,a/2 

M QM {t)= ^* QM (x,t)^ QM (x,t)dx , (31) 

J -a/2 

where ^QM(x,t) is the result obtained from conventional propagation calcu- 
lations, when by conventional we mean the solution of the Hermitian TDSE. 
This enables us to define an effective decay rate even in Hermitian QM based 
on EqE2 

k QM (t) =-j t hxM QM {t) . (32) 

The question we address ourselves with now is how to define the interaction 
region ? or similarly what is a in Eg 1311 ? Is the initial assumption that we 
can define the interaction region as the area between the top of the two barriers 
valid? Here, we are not using the Feshbach formalism and therefore we look for 
a simple universal definition. We propose here to define the interaction region 
(the parameter "a" in our ID case or the vector a in multidimensional case) 
based on the NH calculations. More precisely we wish to define the interaction 
region as the region where the resonances are localized. The continuum states 
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in the same energy range have a very small (almost vanishing ) amplitude in 
this region as can be seen in Fig^ This definition is obviously not an exact one 
and its important to note that using the FP formalism we avoid the need to 
define the interaction region and thus a is only relevant when trying to relate the 
results to conventional QM. The results which will be presented below strongly 
support our conjecture. 

Returning to the problem of the decay of the norm of a WP, we first place a 
gaussian WP as shown in Eg 1261 with a = 3. 87et.it. and fco = 0. It is evident on 
Fig|2]that on the long time scale there is a remarkable correspondence between 
the results of Eq[2HI and Eq|^ On short times there is deviation which results 
from the fact that the decay of the resonances starts at time t — whereas for an 
Hermitian WP it will take time to reach the boundary of the interaction region. 
The idea is simple. Using the NH F-product formalism the resonance states 
decay at any time including the extremely short time regime. In conventional 
propagation calculations the initial WP oscillates in between the barriers and a 
significant part of the tunnelling takes place when the oscillating WP hits the 
inner classical turning points of the potential barriers. Therefore the deviation 
between the results obtained from NH-QM and conventional QM calculations 
is during the time it takes for the initial WP to reach the inner classical turning 
points. To further illustrate this if we study a narrow WP with a = 0.71a.u., 
which is fully localized between the barriers inside the potential in Fig. ^ h 
is obvious from the above argument that in Hermitian QM it will take the 
WP some time to escape out of the barriers while in NH-QM decay starts 
immediately at t = (see FiJ3J- When one studies the behavior of different 
WP's, the best choice for a in EqEH] varies and depends on the WP but it is 
always in the vicinity of the top of the barriers in Fig^ This can be understood 
again by realizing that the continuum wavefunctions under the top of the barrier 
are localized outside of the potential while the resonance states are localized 
inside the potential (even over the barrier), therefore, in NH-QM we have a 
clear distinction between the "system" and its "surrounding". 
In previous nuclear physics studies of the decay of narrow isolated resonances 
which are associated with the NH effective Hamiltonian given in Eq^| the 
decay law has been obtained by the conventional scalar product Namely, 
following this approach the effective decay rate is given by: 

*e//(t) = -|in(**(i)|**(t)) , (33) 

where ^> R (t) is the solution of Eq^ which can be expanded in the basis set of 
the eigenf unctions of EqEl {^a} (with corresponding eigenvalues E a ), such 
that ^/ R (t) = J2 a C R 4> R - Using this approach one gets that, 

M NP = (* R (i)|* R (t)) = "£(C R ,yC R e X p[-<E« - E a ,)t/h] (<f> R \J> R ) , 

a, a 7 

(34) 

where the notation NP stands for the formalism used in nuclear physics in 
Ref.|2Bl- When the resonances are not isolated there is a deviation from the 
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effective decay rate due to interference between the different states. When trying 
to apply this approach to our simple test model, one fails to get converged results 
and moreover the behavior of Mnp based on Eq|M| as can be seen in Fig[3 is 
quite irregular as it oscillates instead of constantly decaying. This shows that 
the formalism depicted in Eq|!2]can only be used in specific cases whereas the 
FP formalism as portrayed in Eq. 1291 should be used in general. 
Let us now test the FP formalism for calculating time dependent expectations 
values for NH Hamiltonians. If one wishes to measure an observable quantity 
A which will be defined in NH-QM, within the FP formalism, by: 

(* FP (t)\A\*FP(t)) J^ FP ( x ,t)A^ FP (x,t)dx 

u (* FP (t)\*Mt)) J^ L FP (x,t)^ P ( x ,t)dx ' (0) 

where in principle A can get complex values. In such a case the phase of A 
should be a measurable quantity (and therefore should be 9 independent al- 
though both ^ FP (t) and ^ FP (t) vary with 8 (the rotational angle associated 
with the CS transformation). In our case (i.e., the Hamiltonian is NH due to 

the CS similarity transformation.) A = SAS^ 1 is the scaled operator for the 
desired observable. As we have seen in Fig's 12131 the NH formalism describes 
the part of the WP which remains in the interaction region. Thus, we expect 
to find correspondence to a quantity similar to that defined in Eq|^ We will 
label < A > as an observable quantity in the interaction region calculated by 
the conventional (Hermitian) quantum-mechanics approach. 

< A(t) >= 'j- , (36) 

where 9(x, t) is the solution obtained by solving the conventional TDSE. Re- 
turning to our one-dimensional test model, the results for the observable A = x 
of the average position for a WP placed at x = with width, a = 3.87a.u., 
and initial momentum k — la.u. based on Ea.'s l3*5l36l are given in Fig0J Once 
again there is very good correspondence on the long time scale while on the 
short time scale there are deviations. This suggests that the NH formalism 
which describes the resonance states applies on times longer than some initial 
rearrangement time in which the resonance states are populated. 
Since the scaling parameter is not associated with any physical quantity, time 
dependent observables should get real values in spite of the analytical continua- 
tion of the Hamiltonian. Indeed, the results presented in Fig|3]show very clearly 
that at any given time the mean position, x(t),{A — x in Ea l35fl . and mean mo- 
mentum p(t) ,(A — p x in Ea l35f) are real quantities even though they have been 
obtained by NH calculation based on Eq|3S] The phase of both the average po- 
sition and momentum is zero with abrupt shifts in 7r every time an observable 
vanishes (and thus can't be measured at that given time). The jumps in tt of 
the phase of A (x(t) or p(t)) result from the fact that the observable changes 
its sign from positive to negative or vise-versa and therefore would appear also 
in standard QM calculations. 
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5 Conclusions 



We have shown that using the recently introduced FP formalism it is possible 
to overcome the time-asymmetry problem in NH-QM and to propagate a WP 
which populates many resonance states. This formalism describes the dynamics 
without the need to separate the entire space into "system" and "surrounding" 
which is necessary when using Hcrmitian QM. The NH description of the system 
is valid on the timescales which are long enough to regard only the localized part 
of the WP after the scattered part of the WP has left the interaction region. 
It has also been shown that observables such as position and momentum in the 
F-product formalism obtain real values despite the NH nature of the Hamilto- 
nian. This reasserts the validity of this formalism which makes it applicable to 
systems where particles are temporarily trapped in the interaction region. The 
application of these novel approach to many electron systems is under current 
study. 
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Figure 1: (Color Online) Eigenfunctions of the model potential in Eq[2Hl The 
solid line is for the bound state; the dot-dashed line is for an isolated resonance; 
The long dashed line is for a continuum state just under the top of the barrier; 
and the dashed line is for a broad resonance over the barrier. 
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Figure 2: (Color Online) The norm of wide Gaussian WP based on the F- 
product formalism as given by Ea l29l (dot-dashed line) as a function of time in 
comparison with the part of the Hermitian WP localized inside the interaction 
region (long-dashed line) as given in EqOH] The dashed line shows the norm in 
the formalism portrayed in Eg 1341 In the inset the initial gaussian WP with 
width of a = 3. 87a. it. and ko = Oa.u. . 
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Figure 3: (Color Online) The same as Fi^21 but with an initial state of a narrow 
gaussian, a — O.Ha.u., fully localized in the interaction region around x = 
(shown in the inset). The dot-dashed line is for Ea 1291 and the long dashed line 
for EqEHl 
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Figure 4: (Color Online) The average position of a gaussian WP with initial 
momentum of k = la.u. and a = 3. 87a. it. as a function of time. The solid 
line is for the solution in the F-product formalism as given in Ea l351 vhile the 
dashed line is for the Hermitian part of the WP localized inside the interaction 
region as depicted in Eq[3f)| The inset shows the difference between the two as 
a function of time. 
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Figure 5: (Color Online) The average position (bottom) and momentum (top) 
of a gaussian WP with initial momentum k = 1 in the F-Product formalism and 
their corresponding phases based on Eq[35] The solid line is for the real part of 
the average position/momentum while their corresponding phases are given by 
the long-dashed lines. The interval of time on which these results appear was 
chosen arbitrarily in order to have clear presentation, but the overall behavior 
is similar in all range of time. Due to the abruptness of the change the shifts in 
pahsesd are not exactly of tt strictly for numerical reasons. 
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